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The simulation of quantum mechanical systems with classical com- 
puters appears to be a computationally intractable problem. In 1982 
FeynmarP reversed this observation, suggesting that quantum mechan- 
ical systems have an information processing capability much greater 
than that of corresponding classical systems, and thus could be used 
to implement a new type of powerful computer. In 1985 Deutsch^ de- 
scribed a quantum mechanical Turing machine, showing that quantum 
computers could indeed be constructed. Since then there has been ex- 
tensive research in this field, but while the theory is fairly well under- 
stood actually building a quantum computer has proved extremely dif- 
ficult, and only two methods have been used to demonstrate quantum 
logic gates: ion traps, and nuclear magnetic resonance (NMR).'^'^ 
NMR quantum computers have previously been used to demonstrate 
quantum algorithms to solve the two bit Deutsch problem.l^l^ Here we 
show how such a computer can be used to implement a fast quantum 
search algorithm initially developed by Grover.^^° 

Among other applications Grover's algorithms enable an extremely rapid 
search over the domain of a binary function to find elements for which this func- 
tion is satisfied (that is, the function has the value 1). This approach is simpler 
if the number of satisfying values is known beforehand, and is particularly simple 
when precisely one quarter of the elements in the domain satisfy the function.^ 
The algorithm can be demonstrated using a computer with two quantum bits 
(qubits) to search a two bit domain in which one of the four elements satisfies 
the function. A classical search of this domain would require between 1 and 3 



evaluations of the function to find the satisfying element, while a quantum search 
can find this element with only one function evaluation. In the more general 
case of searching a domain of size N for one of k satisfying elements the classical 
search requires about ^{N/k) function evaluations, while the quantum searchP^ 
requires only 0{^/N/k). 

A quantum circuit for implementing a quantum search in a two qubit system 
is shown in Fig. ^ This algorithm uses a single function evaluation to label the 
single state which satisfies the function, followed by a series of gates which drive 
the system into this particular state. There are four possible functions / with a 
unique satisfying element, and these functions are conveniently labelled by the 
bit pattern of the satisfying element. The algorithm starts with the quantum 
computer in state |00) and ends with the computer in a state corresponding to 
the bit pattern of the unique element, and so this element can be immediately 
identified by determining the final state of each qubit. 

This algorithm was implemented using our two qubit NMR quantum com- 
puter, described in Ref. This computer uses the two spin states of nuclei 
in a magnetic field as qubits, while radio frequency (RF) fields and spin-spin 
couplings between the nuclei are used to implement quantum logic gates. Our 
molecule, partially deuterated cytosine, contains two nuclei, and thus can be 
used to implement a two qubit computer. The pseudo-Hadamard gates (h) were 
implemented using 90° pulses, while the function evaluation was performed using 
the pulse sequence shown in Fig. |21 The phases of five of the pulses depend on 
which of the four possible functions is to be implemented, and these phases should 
be set as shown in the figure caption. The final gate, Uqq, is easily implemented, 
as it is identical to U . 

The algorithm should start with the computer in state |00), but with an 
NMR quantum computer it is not practical to begin in a true |00) state. Using 
the methods of Cory et al^ it is, however, possible to create an effective pure 
state, which behaves in an equivalent manner. Similarly it is not practical to 
determine the final state directly, but an equivalent measurement can be made 
by exciting the spin system with a further 90° pulse, and observing the phases 
of the resulting NMR signals. The absolute phase of an NMR signal depends in 
a complex manner on a variety of experimental details, and so it is not possible 
to interpret absolute phases, but this can be overcome by measuring a reference 
signal, obtained by applying a 90° pulse directly to the initial state. 

The results of this approach are shown in Fig. IHl Five spectra are shown: a 
reference spectrum acquired using a single 90° pulse, and spectra acquired from 
the same computer implementing the search algorithm for each of the four possible 
functions, /. Each spectrum consists of two closely spaced pairs of lines: each 
pair of lines corresponds to a single qubit, while the barely visible splitting within 
each pair arises from the spin-spin coupling, J, used to implement the two qubit 
gates. To improve the appearance of the spectra the final 90° detection pulse was 
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preceded by a magnetic field gradient pulse, which acts to dephase the majority 
of any error terms which might occur. 

The reference spectrum corresponds to the computer being in state |00), and 
the phase of this spectrum was adjusted so that both lines are in positive ab- 
sorption phase (that is, pointing upwards). The same phase correction was then 
applied to the other four spectra, allowing positive absorption lines to be inter- 
preted as qubits in state |0), while negative absorption lines can be interpreted 
as qubits in state The left hand pair of lines arises from the first spin, and 
thus corresponds to the first qubit, while the right hand pair of lines corresponds 
to the second qubit. Thus, for example, spectrum (c) corresponds to the state 
1 01). Examining the four spectra (b)-(e), it is clear that our implementation of 
a quantum search using the function fab leaves the computer largely in a final 
state \ab), much as expected. 

There are however small but significant errors in the calculation, which result 
in distortions in the final spectra. While most of these distortions are removed 
by the field gradient pulse, their effects remain visible as variations in the heights 
of the NMR lines. These errors arise from a variety of causes, but the most 
significant is errors in the NMR pulse sequence used to implement the calculation. 
These pulse sequences require a large number of selective pulses, that is pulses 
which only affect one of the two spins making up the computer. In practice it is 
difficult to achieve the desired effect at one spin while leaving the other entirely 
unaffected, resulting in errors in the final result. Interestingly the distortions 
are much worse in some cases than in others: they are particularly bad in the 
spectrum obtained when the function is foi- Understanding this variability may 
lead to techniques for reducing the distortion in all cases. 

Implementing Grover's algorithm is a major step forward for NMR quantum 
computing (since this letter was first submitted an implementation of Grovers's 
algorithm using hetero nuclear NMR has been published^, but is by no means 
the limit of what can be achieved. Preliminary studies have been made of systems 
containing three qubits,!^ and it should be possible to build larger NMR quantum 
computers, allowing the implementation of more complex algorithms. 
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Figure 1: A quantum circuit for the implementation of a quantum search al- 
gorithm on a two qubit computer. Gates marked h (pseudo-Hadamard gates) 
act to take a single eigenstate to a uniform superposition of the four possible 
eignestates, while gates marked h"^ implement the inverse operation. The first 
two qubit gate [//^^ corresponds to an evaluation of the function fab, replacing 
an eigenstate \ij) by if i — a and j — b, while C/qo simply replaces |00) by 

-|00). 
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Figure 2: NMR pulse sequence used to implement C//^^. Narrow boxes correspond 
to 90° pulses, while wide boxes are 180° pulses; the upper and lower lines refer to 
the nuclear spins corresponding to the first and second qubits respectively. The 
time period r, during which no pulses are applied, is set equal to 1/4 J, where J 
is the size of the spin-spin couphng between the nuclei. The phase of each pulse 
is written above it, and for five of the pulses this phase depends on which of the 
four functions fab is to be implemented. These phases should be set as follows: 
/oo, = +y, (f) = +x, tjj = -y; /oi, 9 = +y, (f) = -x, i/j = +y; /lo, 9 = -y, 
= -X, -0 = -y; /ii, 9 = -y, (j) = +x, tjj = +y. 
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Figure 3: Experimental spectra from our NMR quantum computer. Spectrum (a) 
is a reference spectrum, used to determine the absolute phases of the NMR signals, 
while spectra (b)-(e) were acquired from the same computer implementing the 
quantum search algorithm using each of the four possible functions: (b) /oo, (c) 
/oi, (d) /lo, (e) fu- These four spectra were processed using the reference phase 
obtained from spectrum (a), and so the phases of the signals can be interpreted 
as states of the corresponding qubits, as described in the main text. 
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